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ABSTRACT 

The power of jets from black holes are expected to depend on both the spin of 
the black hole and the structure of the accretion disk in the region of the last stable 
orbit. We investigate these dependencies using two different physical models for the 
jet power: the classical Blandford-Znajek (BZ) model and a hybrid model developed 
by Meier. In the BZ case, the jets are powered by magnetic fields directly threading 
the spinning black hole while in the hybrid model, the jet energy is extracted from 
both the accretion disk as well as the black hole via magnetic fields anchored to 
the accretion flow inside and outside the hole's ergosphere. The hybrid model takes 
advantage of the strengths of both the Blandford-Payne and BZ mechanisms, while 
avoiding the more controversial features of the latter. We develop these models more 
fully to account for general relativistic effects and to focus on advection-dominated 
accretion flows (ADAF) for which the jet power is expected to be a significant fraction 
of the accreted rest mass energy. 

We apply the models to elliptical galaxies, in order to see if these models can 
explain the observed correlation between the Bondi accretion rates and the total jet 
powers. For typical values of the disk viscosity parameter a ~ 0.04 — 0.3 and mass 
accretion rates consistent with ADAF model expectations, we find that the observed 
correlation requires j > 0.9; i.e., it implies that the black holes are rapidly spinning. 
Our results suggest that the central black holes in the cores of clusters of galaxies 
must be rapidly rotating in order to drive jets powerful enough to heat the intracluster 
medium and quench cooling flows. 

Key words: accretion, accretion disks - black hole physics - galaxies: active - 
galaxies: jets - X-rays: galaxies - MHD 



1 INTRODUCTION 

Recent high resolution Chandra observations of the cores of 
elliptical galaxies provide dramatic illustration of the im- 
pact of jets launched from active galactic nuclei (AGN) on 
the interstellar medium of their host galaxies and/or on 
the intracluster medium if the galaxies reside in clusters. 
These observations have revealed prominent X-ray surface 
brightness depressions corresponding to cavities or bubbles 
inflated by t he jets as they in teract with the surrounding 



2004; [Fabian et all 2006: 



fe.g..lAllen et al.ll2006j. hereafter A06: IB rzan et al 



Taylor et ahlbOOfj ; [Raffcrtv et al 



20061 ). This interaction is presumed to depos it large amounts 
of energy in their environments (e.g., IChurazov et al.l 



E-mail: rodrigo.nemmen@ufrgs.br 



20021: iDalla Vecchia et all 12004 ISiiacki fc Springell 120061 : 
Nusser. Silk fc Babull I 20061) . modifying the hierarchy of 
galaxy formation (e.g., Bower et al,ll2006h and altering the 



evolution of the intracluster medium by counteracting radia- 
tive losses (e.g., McCarthy ct al. 1 12004) . Based on the ener- 
getics of creating the observed cavities, the minimum energy 
associated with the jets range from ~ 10 ergs in galaxies, 
groups, and poor clusters to ~ 10 60 ergs in rich clusters, 
ranking the outbursts as among the most powerful phenom- 
ena in the Universe 

While a detailed understanding of the extragalac- 
tic AGN jet phenomena remains elusive, the combi- 
nation of physical ly insig ht ful analytic studie s (e.g. . 



Blandford & Znajck] Il977i: Blandford fc Pavnel [l982l: 



Begelman et al.lll984l : |Punslv fc Coroni ti 1990; F errarilll998l : 



Meier! I1999T 1200 ll ) and sophisticated general relativistic 
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magnetohydrodynamic (MHD ) numerical simulations (e.g. 
Koide et al.ll2000l; |Koidell2003t iMcKinney fc Gammid 12004 



~ ie y. 

De Villiers et all 120051 : iKomissarovl 120051 : lHawlev fc Krolikl 
20061 ) are beginning to yield important insights. There is 



now a general consensus that jets are fundamentally MHD 
events. 

The currently favoured models presuppose an accretion 
flow threaded by large-scale magnetic fields flowing onto a 
supermassive black hole. In the neighborhood of the black 
hole, the flow settles into a disk-like structure, in which the 
rate of inward flow of matter depends on the efficiency with 
which its angular momentum can be transferred outward. 
This transfer of angular momentum is mediated by MHD 
turbulence. Both MHD turbulence and differential rotation 
of the plasma in the body of the disk generate an d inten- 
sify toroidal mag netic fields l|Balbus fc Hawley||l998h . When 
the pressure associated with the toroidal fields grows strong 
enough, the field lines escape from the disk forming a rotat- 
ing helical tower of field lines above and below the disk. Cen- 
trifugal forces associated with this rotating magnetic field 
helix drives any plasma trapped onto the field lines upward 
and out of the disk, ge nerating outflows. This is the crux of 
the Blandford-Payne l|Brandford fc Pavnel [l982j ) model for 
jets, wherein the magnetic fields extract energy from the 
rotation of accretion disk itself to power the outflows. De- 
pending on the detailed structure of the magnetic fields and 
the accretion disk, this mechanism is expected to generate 
outflows ranging from broad, uncollimated winds to highly 
collimated jets. 

In the event that the black hole itself is spinning, the 
magnetic field can also extract the rotational energy of the 
central black hole to power the outflows. In the Blandford- 



Znajek (BZ) model i|Blandford fc Znaiekl Il977f ), the mag- 
netic fields are assumed to be connected directly to the hori- 
zon of a spinning black hole while in the lPunslv fc Coronitil 
l| 19901 ) model, the magnetic fields associated with jet produc- 
tion are anchored to the inflowing plasma inside the black 
hole's ergosphere. In both these models, the dragging of in- 
ertial frames, relative to an observer at infinity, within the 
ergosphere of the rotating black hole results in a rotating, 
tightly wound vertical to wer of fi eld li nes, and hence pow- 
erful outflows. Moreover, I Meier] fl999) has shown that in 
the Punsly-Coroniti-like models, the differential dragging of 
the frames will also act as a dynamo to amplify the mag- 
netic field at the expense of the black hole's rotational en- 
ergy and that this will have the effect of further enhancing 
the jet power. The results of numerical simulation studies 
ilKoidell2003l; IMcKinney fc Gammiel |2004| ; I De Villiers et al.l 
120051 ; lHawlev fc Krolikl 120061 ) are, in a broad brush sense, 
consistent with the BZ and Punsly-Coroniti-Meier type of 
models, indicating that the power of jets depends on both 
the mass accretion rate as well as the spin of the black hole. 
In fact, iMeied (|200lf ) argues that understanding the 



radio-loud and radio-quiet dichotomy in the QSO popula- 
tion requires the jet power to exhibit such dependencies. 
And while this assertion appears to be indicated by ob- 
servations of black hole c andidate systems in our Galaxy 
i|Cui. Zhang fc Chenl Il998h . very little is known about the 
prevailing conditions underlying the jet phenomena in extra- 
galactic AGNs. The recent Chandra observations of bubbles, 
however, offer a unique opportunity to remedy this. The X- 
ray observations not only provide an estimate of the power 



of the jets emitted by the AGNs but also the rate at which 
matter is accr eting onto the bla ck hole-accretion flow sys- 
tem. Recently lAllen et al.l (|2006l ) (hereafter A06) analyzed 
the data for nine nearby, X-ray luminous elliptical galax- 
ies and found a remarkably tight correlation between the 
Bondi accretion rates and the jet powers. Within the con- 
text of viable models for the jet-black hole-accretion flow 
system, such a correlation not only provides insight into the 
efficiency with which the rest energy of the material accret- 
ing onto the black hole is converted into jet energy but also 
the spin distribution of the black holes powering AGN. In 
this paper, we seek to shed light on these issues. 

We begin by considering two physical models for the jet 
power: the BZ model as d escribed a bove, and a hybrid mode l 
proposed bvlMeierl (|200lT ) (see also. lPunslv fc Coronitilll990l ; 
iMeierl 119991 ) . in which the magnetic field threads the plasma 
throughout — in the body of the disk where frame-dragging 
is negligible as well as within the ergosphere. This lat- 
ter model combines the Blandford-Payne-like disk accel- 
eration mechanism with the Blandford-Znajek-like scheme 
that draws upon the rotational energy of the black hole. 
We improve upon these models by incorporating important 
general relativistic effects that previously we re introduce d 
via highly restricted approximation. Following lMeierji|200ll ). 
we further couple our models for the jet production to the 
ad vection-domina t ed accretion flow model (hereafter ADAF 
B iNaravanl 120051; iNaravan. Mahadevan fc Quataertl Il998l . 
hereafter N98; lNemmen et al.ll2006h . There is a considerable 
body of work indicating that the launching of jets is most ef- 
ficient when the accretion flow is advection-domin ated (e.g., 
iRees et al.lll982l ; lMeierll200ll ; IChurazov et al.lliooBh and that 
jet production is suppressed in the standard thin accre- 
tion disk normally associated with radiatively efficient AGNs 
(|Livio et al.lll999l ; lMeierll200ll ; iMaccarone et al.ll2003T ). 

This paper is organized as follows: In fj2j we describe 
our models for the accretion flow and the jet power. We 
summarize the findings of A06 in Sj3] In we compare the 
correlations between the jet power and accretion rate im- 
plied by our models with the results of A06 and consider 
the resulting implications for jet efficiencies and the dis- 
tribution of black hole spins. We also compare results for 
accretion models with and without winds and explore the 
implications of uncertainties in A06 estimates of jet power 
and mass accretion rates on our findings. Lastly, we highlight 
our main results in [5] and offer some concluding remarks. 



2 MODELS OF ACCRETION FLOW AND JET 
POWER 

2.1 Accretion flow structure 

There are several lines of evidence that suggest that the 
accretion flow onto AGNs with strong jets is best de- 
scribed as an ADAF: First, the observed bolometric (radia- 
tive) luminosities of radio-loud AGNs are typically many 



1 Improved models of advection-dominated accretion flows incor- 
porating winds and convection are also more generally referred to 
as radiatively inefficient accretion flows (RIAF). We simply use 
the original acronym ADAF. 
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orders of magnitude smaller than the luminosities ex- 
pected if the mass were flowing onto t he central black 
holes along thin accre tion disks (e.g., iDi Matteo et al.l 
120031 ; iTavlor et all [2006 ) ; in the latter case, the bolomet- 
ric luminosities correspond to roughl y 10% of the rest- 
mass energy of the accreting matter (Sh akura fc Sunvaevl 
1 19731 ). ADAFs can not only account for the observed 
radiative quiescence of the jet-emitting AGNs but also 
reproduce their nuclear X-ray luminosities with accre- 
tion disk mass flow rate cl comparable to their Bondi 

l2003l ; lLoewenstein et al.ll200ll ; 



rates (|Di Matteo et 
|Pellegrinill2005l ; ISoria 




2006h . The exceptional systems 



(such as the radio-loud quasars) that are both highly lu- 
minous and radio-loud are currently an intriguing puzzle: 
possibly they are high accretion rate systems in which the 
disk is puffed up by t he trapped energy despite the sys- 
tem's low opacity (e.g.,|M cicr 200 11; iMaccarone et alj 2003; 



Kording. Jester fc Fenderl 120061 ; ISikora. Stawarz fc Lasotal 
20071 ). 



Second, the magnitude and the structure of the mag- 
netic fields associated with ADAFs are much more con- 
ductive to the extraction of spin energy from the hole 
than those associated with standard thin disks (e.g., 



Rees et al.l Il982l ; iLivio et alj 1 19991 ; lArmitage fc NataraiaiJ 
19991 ). Third, X-ray binaries with black hole candidates 
(XRBs) in the "low/hard" state display a strong corre- 
lation between the X-ray and radio emissions, which can 
best be understood in the context of an ADAF-jet system 
i|Gallo et al.ll2003l ; iMerloni et ai1l2003l ; iFalcke et al.ll2004 ). 
Low-luminosity AGNs with jets show similar correlations 
and hence , are thought to be scaled-up versions o f the galac- 
tic XRBs ((Merloni et al.ll2003l ; lFalcke et al.ll2004 ). One char- 
acteristic of the ADAF-jet system is low mass flow rates. The 
observations suggest that at high flow rates, the black hole 
systems switch to the "high/soft" stat e and the jet activit y 
is suppressed (|Maccarone et al.l 120031 : [G rccne et alj 2006). 
This observed "quenching" of the jet (as well as the tran- 
sition from the "low/hard" to the "high/soft") is believed 
to be due to a change from an ADAF (powerful jets) at low 
mass flow rates to a thin accr etion disk st ructure (weak jets) 
at high mass flow rates (e.g., Meier 2001). 

We describe th e structure of the AD AF using the self- 
similar equations o f lNaravan fc~Yil ij 19951 ). We are especially 
interested in the analytical equations (see Appendix I) that 
describe the vertical half-thickness of the disk H , the angular 
velocity of the disk Q.' and the magnetic field strength B 
near the black hole in terms of radius R, black hole mass 
M., accretion rate onto the black hole M and advection 
parameter / (assumed « 1). We require these as inputs to 
our models for the jet powe r described in i ]2.2l and l2.3l below. 
We note that although the iNaravan fc Yin 19951 ) equations 
are based on an early model for the ADAF that did not allow 
for mass loss in the form of winds from the accretion flow, 
we have verified that these self-similar solutions provide a 
good approximation to the structure of the inner regions of 



2 We use the term 'mass flow rate' to distinguish the mass flowing 
through the accretion disk from the mass that is accreted onto the 
black hole itself. Due to the mass/energy carried away by winds 
and jets, these two rates may differ significantly. 



ADAF models, like that of iBlandford fc Begelmanj (|l999l) 
(hereafter BB99), that do allow for such losses. 

The equations for H, B and f2' also depend on the 
properties of the accreting fluid, such as its adiabatic in- 
dex 7, its viscosity parameter a and the ratio of gas to 
magnetic pressure /3 in the fluid. These quantities are not 
independent. The value of 7 depe nds on /3 via th e re- 
lationship 7 = (5/3 + 8)/3(2 + /3) (|Esin et all Il997l ) and 
based on the MHD numerical simulations of the evolu- 
tion of the magnetorotational inst ability in accretion disks 
l|Hawlev. Gammie fc Balbuslll995l ). a » 0.55/(1 + /3). We 
assume that the magnetic pressure is related to the field 
strength as P mag = B 2 /8tv. 

We constrain the value of a from recent MHD numerical 
simulations of radiatively inefficient accretion flows, which 
take into account self-consistently the role of the Maxwell 
stresses in establ ishing its value. The general relativistic 
simu l ations (e.g., McKinnev fc Gammiej |2004| ; iHirose et al.l 
l2004lHawlev fc Krolikll2006l ) find that in the inner portions 
of the disk /3 ~ 1 — 10. Although a is not explicitly quoted, 
we can estimate a from these values of /3 using the relation- 
ship between the two, obtaining a ~ 0.04 — 0.3. For com- 
putational purposes, we will derive results for a — 0.04 and 
0.3. To put these a-values in context, w e note that the sim- 
ulatio ns around Schwarzschild holes of iProga fc Begelmanl 
|2003l ) suggest that near the innermost stable circular or- 
bit, q reaches high values: a ~ 0.1 — 0.7. Moreover, recent 
ADAF models of XRBs require values of a » 0.25 in orde r 
to account for the observations (|Quataert fc Naravanlll99^ ) . 

To the above model for the ADAF, we introduce three 
important modifications to take into account general rela- 
tivistic effects induced by the Kerr metric: First, we take as 
input for our jet model the values of H, B and Q' evalu- 
ated at Rm S , the radius of the marginally stable orbit of the 
accretion disk. This radius depends sensitively on the dimen- 
sionless black hole spin parameter j = J/J max ("o/M" in ge- 
ometrized units, a is the specific angular momentum), where 
J is the angular momentum of the ho le and J ma x = GMj / c 
is the maximal angular momentum l|Bardeen et al.ll"l972l ). 
Second, an observer at infinity will see the disk and the 
magnetic fields near the black hole rotate, in the Boyer- 
Lindquist coordinate system, not with an angular velocity 
il' but fl = O' + oj, where uj = —g,f>t/g<p4> is the angu- 
lar velocity, in the same coordinate system, corresponding 
to the loca l spacetime rotation enforced by the spinning 
black hole (|Bardeen et al.l fl972T ) . And third, we take into 
account the field-enhancing shear caused by frame-dragging 
when calculating the magnetic field stre ngth in the inner 
region of the disk, as first suggested by iMeierl (|l999l ) and 
tentatively observed in recent MH D num e rical simulations 
l|Hawlev fc KrolfldbOOrj) . Following IMeierl (|200ll ). we relate 
the amplified, azimuthal component of the magnetic field 
to the unamplified magnetic field strength derived from the 
self-similar ADAF solution as B$ = gB, where g = Sl/fl' 
is the field-enhancing factor. The amplitude of this factor 
depends on the black hole spin through the angular velocity 
of spacetime rotation, uj. In the case of a non-rotating black 
hole, uj — and g = 1 (i.e., no field enhancement). 

The potential imp ortance of th ese effects has been rec- 
ognized previously and|Mcicr (2001) even incorporated them 
into his jet power model, albeit in the form of simplifications 
that blunted their impact. For example, he took g to be a free 
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parameter rather than explicitly relate it to spacetime rota- 
tion induced by the rotating black hole, he adopted a simple 
approximation for ui, and he evaluated his jet model only for 
two limiting radii: 7GM,/c 2 (corresponding to j ~ 0) and 
1.5GM./c 2 (j » 1). In our approach, we are able to capture 
more fully the dependence of the important model parame- 
ters on j and by doing so, are able to explore the sensitivity 
of the accretion flow and the jet power solutions to variations 
in the value of the bla ck hole spin param eter. 

Finally, we follow iLivio et al.l (|l999T l in assuming that 
the poloidal and azimuthal components of the magnetic field 
are related to each other as B v ~ H/R B^, where R is the 
radius in cylindrical coordinates, a relationship based on the 
assertion that the strength of the poloidal field is limited by 
the vertical extent of turbulent eddies in the disk. In the 
case of an ADAF, this yields B p as B^ because H ~ 7?. 



2.2 The Blandford-Znajek jet model 

According to the BZ model, the total power of the re- 
sulting jet is give n by (e.g., iMacdonald fc Thornd Il982l : 
iThorne et al]|l986l ): 



r,BZ 1 2 D 2 n 2 .2 
Pjet = -^^fBj_R h J C, 



(1) 



where Rh = [1 + (1 — j 2 ) 1 ^ 2 ]GA/./c 2 is the horizon radius, 
B± is the strength of the magnetic field normal to the hori- 
zon, and the factor u>f = Uf(fiH — Q,f)/&h depends on 
the angular velocity of the field lines Of relative to that 
of the hole, Qh- Following the usual practice, we assume 
that ujf = 1/2, which maximizes the power output (e.g., 
Macdonald fc Thorne|[l982l ; IThorne et alJll986h . As for B ± , 



Livio et all l| 19991 ) has argued that the field threading the 



horizon ought to be comparable in strength to the field 
threading the inner regions of the accretion flow. Hence, 
we take the field strength at the horizon to be the same 
as that at the radius of the marginally stable orbit of the 
accretion disk; that is, B± sa B p (R nls ) g(R lnB )B(R lnB ). 
These field configurations are consistent with those seen 
in the in the jet launch regions in nu merical simulations 
of accretion onto Kerr blac k holes (e.g., iHirose et al. | |2004 
iMcKinnev fc Gammidl20oll ). Finally since the accretion rate 
enters the jet power only through the field strength, the 
jet power only depends on the accretion rate measured at 
the marginally stable orbit of the accretion flow Af ms = 
M(R ms ). If the mass is cons erved in the disk as i n the early 
version of the ADAF model (|Naravan fc Yill 19951 . N98) then 
M is constant with radius, but contemporary ADAF models 
often allow for mass loss from the accretion disk in the form 
of winds and outflows (see £|2.4[1 . 

The primary advantage of the BZ model is its simplic- 
ity. The picture describes a black hole whose event hori- 
zon is equivalent to a rotating conducting surface with sur- 
face charges and currents. In effect, the horizon is a unipo- 
lar inductor onto which the magnetic fields are attached. 
The magnetic fields are torqued by the rotation of the black 
hole and in turn, drive an outflow. There ar e , how ever, con- 
cerns about the model. IPunslv fc Coronitil (Il990l ) have ar- 
gued that the very notion of flux being emitted from the 
event horizon renders the model unphysical because the 
event horizon is casually disconnected from events upstream. 



Disturbances generated at the event horizon can only prop- 
agate inward. Moreover, the BZ model is also formulated 
on the assumption of zero accretion that underlies the BZ 
model; the BZ model is deriv ed in the limit of vani shing 
plasma density and pressure l|Hawlev fc Kr olik 2006). On 
the other hand. IMcKinnev fc Gammiel 1 20041 ) have reported 
good agreement between several aspects of the BZ model 
and their numerical simulation results, and for this reason, 
we consider t his model in the present wo rk. We do note, 
however, that IMcKinnev fc Gammiel (|2004l ) also reported a 
growing difference between the model and their results with 
increasing black hole spins, a result that can partly be at- 
tributed to the fact that the BZ model is derived in the limit 
of slowly rotating black holes. 

2.3 Hybrid model 

Contrary to the premise of very low plasma densities that 
underlies the BZ model, accretion of matter is presumed 
to be a key element of real AGN systems. In fact, numer- 
ical simulations show that the coupling between the accre- 
tion flow and the magnetic fields is an essential element of 
jet production. For this re ason, we cons ider a second jet 
model, the hybrid model of iMeierl (|200ll ). As noted previ- 
ously, this model is constructed so that large-scale magnetic 
fields thread the accretion disk outside the ergosphere as 
well as rotating plasma within the ergosphere that is flow- 
ing onto the black hole. Hence, the model is able to draw 
upon both the rotational energy of the accretion disk as well 
as the spinning black hole in order to drive outflows, though 
the extraction of the black hole rotational energy occurs in- 
directly through field lines anchored to the plasma subject 
to spacetime rotation imposed by the black hole. This model 
also takes into account the effects of field amplification by 
both the differential rotation of the plasma in the body of 
the disk and the differential frame-dragging. As discussed by 
IMeierl (|200ll ). the jet power in this model is a strong func- 
tion of the thickness of the accretion disk and the black hole 
spin; strong magnetic fields and rapid rotation, the neces- 
sary ingredients for the launching of powerful jets, only arise 
when the disk is thick and the hole is spinning rapidly. 

Following IMeierl l|200ll ). the total jet power for the hy- 
brid model is given by 



pdisk 
-fjct 



(B^HRQ) /32c, 



(2) 



where B^ — gB and Q = SI'+oj. All quantities are evaluated 
at R = ii m s, which is also assumed to be the approximate 
characteristic size of the jet-formation region. We note that 
we denote the jet power for this model with the superscript 
"disk" in order to highlight that in addition to drawing upon 
the rotational energy of the spinning black hole, this model 
also draws energy from the accretion disk. 

As a matter of interest, we note that there is consid- 
erable support for the Punsly-Coroniti-Meier component of 
t he hybrid model from recent numerical simulations studies 
llKoidell2003l ; IMcKinnev fc Gammiel liooi | De Villiers et al.l 
120051 ). On the other hand, the relevance of the Blandford- 
Payne mechanism, which describes the extraction of jet 
power from the rotation of the accretion disk and which 
is incorporated in our hybri d model, is very much a mat- 
ter of debate. For example, lHawlev fc Krol ik (2006) have 
recently argued that this mechanism is not at all important 
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Figure 1. The jet powers for our modified Blandford-Znajek 
(Equation [TJ and the hydrid model (Equation [2j as a function 
of black hole spin parameter j. The curves are computed assum- 
ing a fiducial value of M ms = 10 -3 Mq yr — 1 , and for two different 
values of the viscosity parameter: a = 0.04 and 0.3. The right axis 



shows the jet efficiency rjj et = 
scale linearly with M ms , Jjjet 



= Pj t/A/ ms c . Since the jet powers 
depends only on j and a. 



for both the BZ and the hybrid models is a strong func- 
tion of the black hole spin, spanning a range of three or- 
ders of magnitude as the black hole spin varies from to 
1. This strong dependence is largely due to our improve- 
ments over previous models, which di d not incorporate care- 
fully the physics of the Kerr metri c (IGhosh fc Abramowicj 
1 19971 ; lArmitage fc Nataraianlll999l ; lMeierll200ll ). This type 
of ji-dependence for the jet power has previously only been 
seen in complex numerical MHD simulati ons of jet formation 
|McKinnevll2005l ; lHawlev fc Krolikll2006h . 

The right axis of Figure [T] shows the jet efficiency fac- 
tor, defined as ?7jet = Pj c t/M ma c 2 . Since our model jet powers 
scale linearly with M ms , the corresponding jet efficiency fac- 
tor depends only on the black hole spin parameter j and 
the viscosity parameter a. For high values of j, the effi- 
ciencies predicted by our two jet models are comparable to 
the Novikov-Thorne thin disk radiative efficiency, which is 
related to t he binding energy of the innermost stable cir- 
cular orbit (|Novikov fc Thorndll973l ). and the jet carries 
away a significant fraction of the rest mass energy flowing 
through the disc. We note that for the max imal black hol e 
spin, j — 0.998, as implied by the analysis of lThorriel (| 19T4T ) . 
the jet efficiencies are rj = 0.22 (BZ model) and 0.48 (hybrid 
model) for a = 0.04. On the other hand, if a = 0.3, the 
efficiencies drop to r\ — 0.07 (BZ model) and 0.24 (hybrid 
model) . 



for understanding the collim ated jet outflow s seen in their 
numerical simulations while iBlandfordl (|2005h has suggested 
that such conclusions are premature and that the simulation 
results are dependent on the initial conditions adopted. 

2.4 Jet models properties 

There are several interesting properties of our jet models 
that ought to be highlighted. First, the jet power for nei- 
ther the BZ nor the hybrid models (Equations Q] and re- 
spectively) depends on the black hole mass. Moreover, the 
dependence on the accretion rate enters the jet power only 
through the field strength (see Appendix I) and for the com- 
bined accretion flow-jet power model under consideration 
here, the jet powers for both models can be expressed as 
Pj et (a,j, Mms) oc M ms ; that is, the jet power depends lin- 
early on the mass accretion rate. However, both have a com- 
plicated dependence on j and a (see Appendix I). Figure [T] 
illustrates the spin dependence of both the BZ and the hy- 
brid Pj c t models, for viscosity parameters a = 0.04 and 0.3, 
and a fiducial value of M ms = 10 -3 Mq yr -1 . This corre- 
sponds to one tenth of the mean Bondi accretion rate of the 
elliptical galaxies studied by A06. 

As Figure [T] shows, both the BZ and the hybrid Pj e t 
models behave similarly and give comparable results for in- 
termediate and high values of j. The slight difference in the 
jet power between the two models is not significant given 
the uncertainty in our understanding of the detailed mecha- 
nisms underlying jet formation. At low spin values, however, 
the jet power in the hybrid model is considerably higher 
than that for the BZ case. The drop-off in the BZ case re- 
flects the decline in the black hole rotational energy avail- 
able to power the jet while in the hybrid case, the decline is 
limited by contributions from the accretion disk. The most 
important point made by Figure [1] is that the jet power 



3 THE EMPIRICAL M-Pjet CORRELATION 

Given a distribution of gas about a central black hole, the 
most simple configuration describing the accretion of the 
gas o nto the the black hole is the Bondi flow model (|Bondil 
ll952T ). which assumes a non-luminous central source and a 
spherically symmetric flow with negligible angular momen- 
tum. The resulting Bondi accretion rate can be written as 
Afeondi = n\c s pr A , where ta = 2GM,/c 2 s is the accretion 
radius, G is the gravitational constant, M. is the black hole 
mass, c s is the sound speed of the gas at ta, p is the density 
of gas at ta and A is a numerical coefficient that depends on 
the adiabatic index of the gas. This estimate is frequently 
used in studie s of the central X-ray emitting gas in el liptical 
galaxies (e.g.. |Pi Matteo et al.ll2003l ; |Pellegrini||2005l ). 

Recently, A06 derived the Bondi accretion rates and 
jet powers from Chandra X-ray observations of nine nearby, 
X-ray luminous giant elliptical galaxies that show evidence 
of jet-inflated cavities in their central regions, and found 
a tight correlation between the two. The jet powers were 
estimated from the energies and timescales required to in- 
flate cavities observed in the surrounding X-ray emitting gas 
such that Pjet = -E/tagc, where E is the energy required to 
create the observed bubbles and i ag e is the age of the bub- 
ble. A06 estimates are based on the assumption that the 
X-ray bubbles are inflated slowly. The values of MBondi for 
the systems were calculated from the observed gas temper- 
ature and density profiles. In most cases, the bondi radius 
is not observed directly, and the appropriate density and 
temperature are determined by extrapolating the observed 
data, typically by a factor of 3 or greater, in radius (cf., 
iRaffertv et aLll200rJ ). The black hole masses were deduced 
from the optical velocity dispersion measurements using the 
correlation between central black hole mass and velocity dis- 
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implies that a non-negligible fraction (Pj e t/(MBondiC 2 



-1.5 -1.0 -0.5 0.0 0.5 1.0 
log(P jet /1 43 erg s" 1 ) 

Figure 2. The empirical relationship between the Bondi accre- 
tion power (PBondi = O.lMeondiC 2 ) an d jet power (-Pjet) f° r nine 
nearby, X-ray luminous giant elliptical galaxies derived by A06. 
The error bars plotted include a systematic uncertainty of 0.46 
dex in log PBondi- The fitted power-law model predicted by our 
jet models (A = 0.69 ± 0.19 and B = 1) is represented by the 
dashed (best-fit) and dotted (error bars) lines; the best-fit power- 
law model determined by A06 is shown as the solid line. 



persion of lTremaine et all l|2002l ). Since A06 presents a state 
of the art analysis of current data, we accept the estimates of 
the Bondi accretion rate and jet power at face value in what 
follows. Clearly it will be important to see how these obser- 
vational constraints improve with future X-ray missions. In 
section 14.21 we consider the impact of uncertainties in the 
calculation of the jet power from the physical properties of 
cavities observed in the Chandra data. 

The A06 results for the nine systems are shown in Fig- 
ure [2] where PBondi is the total accretion power released 
for an efficiency of 10%, PBondi = 0.1A/ Bo ndiC 2 . The error 
bars plotted include a systematic uncertainty of 0.46 dex in 
log PBondi, implied by the intrinsic dispersion of 0.23 dex in 
the log M, — log a relation (see A06) . 

The correlation between A^Bondi (or equivalently, 
PBondi) and jet power Pj e t was expressed by A06 as a power- 
law model of the form 



log 



Pbo 



10 43 erg s _1 



= A + B log 



10 43 erg s" 



(3) 



with A = 0.65 ± 0.16 and B = 0.77 ± 0.20. This power-law 
model is shown in Figure [5] as a solid line. 



4 Af— Pjet RELATION: IMPLICATIONS FOR 
BLACK HOLE SPIN AND ACCRETION 
RATE 

The A06 results are interesting in two regards: First, they 
show that the central AGNs in these systems are extremely 
sub-Eddington; that is, they have extremely low accretion 
rates compared to the Eddington rate, Maondi/AfEdd < 
10" 3 , where Af Edd = 22A/./(10 9 M©) M©yr _1 . Our ADAF 
model for the accretion flow is only valid for such highly sub- 
Eddington flows. Second, as noted by A06, the correlation 



2.2l L n "% for P ic 



10 



43 erg s 1 " 



of the energy associated 



with the rest mass of the gas entering ta is channeled into 
jet power. 

4.1 Reconsidering the mass accretion rate 

Presumably less matter than the amount predicted by the 
Bondi rate gets down to the black hole for several reasons. 
As the gas in an ADAF has angular momentum, accretion is 
driven not just by gravity as in the Bondi flow but also by the 
rate of angular momentum transport characterized by the 
parameter a in our ADAF model; given a certain ambient 
density in the external medium, the accretion rate predicted 
by the ADAF model is lower tha n the Bondi accretion rate : 
A/adaf ~ aA/Bondi (e.g., N98; iProga fc Begelmanl [20031 ) . 
with a < 1. Furthermore, some part of the gas may be pre- 
vented from being accreted due to winds and/or convection 
(e.g BB99;lQuataert fc Gruzinovll2000l ; [proga fc Begelmanl 
120031 ; llgumenshchev et al.ll2003h occuring in the ADAF, re- 
ducing even more M ms compared to Maondi- We allow M ms 
to be smaller than MB on di by introducing the parameter 
fBondi such that M ms — ^Bond iAfaondi- The parameter ^Bondi 
represents the fraction of material supplied by the external 
medium at the Bondi radius that ultimately reaches the in- 
nermost stable circular orbit of the accretion flow and gets 
accreted afterwards. Therefore this parameter encompasses 
our ignorance about the possible physical processes that may 
modify the density profile of the accretion flow with respect 
to the ADAF solution. We note that here we are only con- 
sidering the simplest situation where the density profile is 
affected. Winds too will reduce the rate of energy and an- 
gular momentum accretion onto the central black hole. We 
defer the discussion of such complications to H4.3I 

Including this modification into our models gives jet 
power relations of the form 



Pjct ondi? ^Bondi? jf) ^ -^Bondi- 



(4) 



This relationship has the same functional form as the corre- 
lation found by A06 (Equation |3J) and we can now use the 
above equation to constrain the main model parameters (a, 
eBondi and j) using A06 results. 

Comparing Equations [3] and U it follows that our mod- 
els predicts the slope B = 1, which is somewhat higher than 
the value obtained by A06, B = 0.77 ± 0.20, although the 
difference is not statistically significant (see Fig. [2}. As A 
and B in Equation[3]presumably are strongly correlated, we 
fit A06's data to a power-law model with the fixed value 
B = 1 to find the corresponding value of A. Using the x 2 
fit statistics, which accounts only for errors in the values of 
PBondi, we find A = 0.69 ± 0.19 with X 2 = 2.4 for 8 degrees 
of freedom, which indicates that our B = 1 fit also provides 
a good description of the data, and hints that the system- 
atic errors in the data points are not fully independent. The 
power-law model with A — 0.69 and B = 1 is plotted in 
Figure [2] as the dashed line; the dotted lines delimitate the 
corresponding statistical uncertainty in this power-law. 

In the B = 1 model, the value of the parameter A is a 
measure of the jet efficiency r/j a t. Setting B = 1 in Eq. [3] we 
have 



( = [ ^ondi | _ 



V Pi 



jet 



0.1 



^jotEBondi 



(5) 
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Figure 3. Relation between the parameters r)j et and eBondi re- 
quired to reproduce the observed correlation between -MrSondi an d 
Pjet- Each line correspond to a different value of A (shown above 
the lines), the thicker line represents the best- fit value of A. We 
limit the value of efficiency to the theoretical maximum for the 
hybrid model with a = 0.04 (48%). The dotted line corresponds 
to the maximum efficiency for the Blandford-Znajek model with 
a = 0.04. These correspond to maximally spinning black holes 
with j = 0.998. The inferred efficiencies decline with increasing 
a since the plausible values of eBondi increase. For a = 0.3, the 
maximal efficiencies are r] = 0.07 (BZ model) and 0.24 (hybrid 
model). 

which can be recast as 

log(»?jct) = -(1 + A + loge B ondi). (6) 

Figure [3] shows the jet efficiency factor corresponding 
to the observed value of A = 0.69 ± 0.19 (B = 1), as a 
function of the mass accretion parameter eBondi- The value 
of r/jet » 2% derived by A06 is, in fact, the minimum value 
and corresponds to eBondi = 1 or M ms = Msondi • Since eBondi 
is expected to be significantly less than unity, the required 
jet efficiency will be higher. 

Since in our models for the jet power, the jet efficiency 
factor is a function of the viscosity parameter a and the 
black hole spin j, the above equations imply that the value 
of the parameter A is a function of the parameters a, eBondi 
and j. And therefore, for a given value of a, we can use the 
observed value of A = 0.69 ± 0.19 to explore the range of 
possible values of eBondi and j. As described previously, var- 
ious different analysis and arguments suggest that a ranges 
between 0.04 and 0.3. In the discussion below, we will derive 
results for the two values that bound the range: a = 0.04 
and 0.3. 

Figure [4] shows the parameter space (j, eBondi) corre- 
sponding to the observed values of A for the BZ (dashed 
line, jgUl and the hybrid (solid line, ffO]) jet models. The 
two panels show the results for the two different values of 
a near the innermost stable circular orbit, a = 0.3 (/3 ~ 1, 
left panel) and a = 0.04 (/3 ~ 10, right panel). We remind 
the reader that j3 is the ratio of the gas pressure to the mag- 
netic pressure near the radius of the marginally stable orbit. 



The label beside each contour is the value of A for that line. 
We show the effect of considering an uncertainty of ±0.19 
in A. Note that the range we have adopted assumes that 
the uncertainty in each data point is independent. A would 
be smaller if, for example, the jet power was systematically 
under-estimated in all systems (see N4.2|) . 

Based on the arrangement of the contours in this figure, 
we find that the observed tight correlation between accre- 
tion rates and jet powers implies a narrow range of black 
hole spins for the elliptical galaxies of the sample of A06, ir- 
respective of the value of a and the specific model adopted, 
with the main result being that the central black holes pow- 
ering the jets in these systems must be rapidly spinning. 
Specifically, we find that if a ~ 0.3 (Figure^), as required 
to ensure an agreement betwe en ADAF model prediction s 
for XRBs and the observations |Quataert fc Naravan|[l999l ) , 
then the hybrid model requires j > 0.75, while the BZ model 
implies j > 0.88. If the value of the viscosity parameter is on 
the low side, a « 0.04 (Figure |4p), both jet models require 
j > 0.75. In all cases, this corresponds to rjjet > 2%. 

These lower limits for both j and rjj ct correspond to the 
situation where rate of mass accretion onto the black hole 
equals the Bondi rate (esondi ~ 1). As discussed above, only 
a fraction of the Bondi rate is likely to make its way to the 
black hole and realistically, M ms < Madaf- This restricts 
the allowed region of parameter space in Figure [4] to that 
below the horizontal dotted lines (eBondi = a) corresponding 
to Mm B = Madaf- 

Subject to these constraints on the mass accretion rate, 
the observed correlations between mass accretion rate and 
jet power imply that in the case of the hybrid model, the cen- 
tral black holes must have spins j > 0.89 if a ~ 0.3. The cor- 
responding range for the mass accretion rate onto the black 
hole is M ms w (0.05 - 0.3)M Bo ndi ~ (0.16 - 1.0)A/ A daf. 
If a is on the low side, then the black holes in the hy- 
brid model must be spinning rapidly (j > 0.98) and the 
corresponding range for the mass accretion rate is A/ ms ~ 
(0.02-0.04)M Bon di « (0.5-1.0)Madaf- For a given value of 
j, the corresponding jet efficiency factor can be read off from 
Figure [3] For example, if a = 0.3 and j = 0.89, ?7j c t = 4.5%, 
giving an overall efficiency of JfteteBondi = 1.3%. 

In the case of the BZ model, the spin distribution cor- 
responds to j > 0.98 (a = 0.3) and j > 1 (a = 0.04). The 
latter is unphysical and indicates that values of viscosity 
parameters as low as a = 0.04 are disfavored. For this line 
of reasoning, we find that the BZ model requires a ^ 0.2 
in order to guarantee j ^ 1 and M ma ^ Madaf- It is also 
possible these "problems" reflect the limitations of the BZ 
model. As indicated previously, the BZ model is, strictly 
speaking, applicable only in the limit of slow rotation where 
the unipolar approximation is reasonable. For this reason, 
and for the clarity of the discussion, we shall hereafter only 
discuss the results in the context of the hybrid model though 
we shall continue to show the results for both models. 

Since our results indicate high spin values for jet- 
emitting black holes, it useful to bear in mind that re- 
ce nt numerical re l ativis tic MHD simulations of thick-disks 
of iGammie et all l|2004h suggest that the accreting plasma 
will bring the black hol e in spin equil ibrium not at the maxi- 
mal value of j = 0.998 (|Thornelll974h but at a lower value of 
j ~ 0.93. This is not a trivial difference. At high spin rates, a 
7% reduction in j translates into a factor of ~ 3 reduction in 
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Figure 4. Range of values for the parameters eBondi and j from the jet models, which reproduce the measured values of A of the 
empirical correlation between MB on di and Pj et (Equation [5j . The dashed lines show the predictions of the Blandford-Znajek model and 
the solid lines show those of the hybrid model. The left panel corresponds to a = 0.3 and the right panel to a = 0.04, the two values that 
bound the range of plausible values for the viscosity parameter. The label beside each line is the value of A for that line. The thickest 
lines correspond to the best fitting value of A. The horizontal dotted lines indicate the value ^Bondi = a (i- e -i W ms = Madaf)- 



the jet power. In the specific case of the A06 elliptical galax- 
ies, we find that black holes in these galaxies must be fed 
at rates tBondi > 0.05, depending on the value of a and the 
mechanism of jet powering at work. In the particular case of 
a hybrid model with a ~ 0.3, the black holes cannot be fed 
at rates much lower than Madaf (i.e, 0.2 < eBondi < 0.3, 
where the lower limit is due to the reduced efficiency asso- 
ciated with j = 0.93 and corresponds to M ms = Madaf), 
limiting the potential role of mass loss through winds. 



4.2 Reconsidering the jet power estimates 

A06 uses the X-ray observations of the intracluster medium 
cavities to estimate the jet power. To derive the energy E re- 
quired to create the observed bubbles in the X-ray emitting 
gas, A06 assume that the cavities are inflated slowly and 
obtain the relationship E — 4PV for 7 = 4/3 (relativistic 
plasma), where P is the thermal pressure of the surround- 
ing X-ray emitting gas, V is the volume of the cavity and 
7 is the mean adiabatic index of the fluid within the cavity. 
This is the minimum energy required to inflate the cavi- 
ties and does not take into account any additional energy 
that may have gone into heating the intracluster medium. 
A more realistic scenario is likely to involve overpressurized 
bubbles, which upon injection expand rapidly to reach pres- 
sure equilibrium with their surroundings, and in the process 
gen erate weak shock wav e s tha t heat the X-ray emitting 
gas. iNusser. Silk fc Babull (120061 ) calculate the total energy 
E' deposited by the jets in such cases and find that 



E' « 3PV 



1/4 



(7) 



gests that if a bubble is injected with overpressure Pi / P > 
10 then E' > 6PV . We consider the possibility that as an ex- 
treme case the bubble energies are twice the value assumed 
by A06 (E' — 8PV), and calculate the impact of this on 
the observed jet powers and on the results derived from our 
models. 



Following A06, 



calculate the modified jet powers 



where Pi is the thermal pressure of the surrounding gas at 
the location where the bubble is injected. Equation [7] sug- 



as ^j'et = E' /t^gc = 2Pj Gt and refit the power-law model 
with B = 1 to the modified data, taking into account the 
proper error propagation on the values of Pj ct - We obtain 
A = 0.39 ± 0.19 with X 2 = 2.4 for 8 degrees of freedom. 
As the assumed value of E' implies higher jet powers, our 
models need higher values of j to reproduce the increased 
values of -P/ e t- In particular, if a ~ 0.3 these values of A im- 
ply jmin ~ 0.84 for the hybrid model and j m - m ~ 0.95 for the 
BZ model. Therefore, a narrower range of still large spins is 
required to explain PL t and the implied numerical values of 
eBondi and cadaf also increase slightly. 

One important uncertainty in the above calculation 
(and also in that of A06) is the estimate of the age of the 
bubbles. A06 calculate the ages using the formula t Cs = 
D/c a , where D is the distance of the bubble centre from the 
black hole and c s is th e adia batic sou nd speed, but as dis - 
cussed bv lBirzan et al.1 l|2004h (see also lRaffertv et al.ll2006h 
there are two other ways o f estimating the ag e of the cavi- 
ties (Equations 3 and 4 of iBirzan et al.l [20041 ) which result 
in longer timescales when compared to t Cs . An underesti- 
mate of the timescales involved impacts the estimate of jet 
power, as determined above, by artificially enhancing the 
jet power above the "true" value. To take this into account, 
we consider the p ossibility that the A06 ages are too low 
by a factor of - 2 (|Birzan et al.ll2004 iRaffertv et"al1l2006l ). 
such that the modified age is t agc « 2t Cs , implying a smaller 
jet power Pj" t = l/2Pj e t- With reference to the original jet 
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power estimates of A06, we find that revising the timescales 
as suggested here results, not surprisingly, in a systematic 
decrease in lower limits of ji-values. For example, in the case 
of the q = 0.3 hybrid model, we find that j m i n « 0.6. 

In spite of the above-mentioned concerns with the 
A06 estimates of the total jet energy and the associated 
timescales, it seems likely that their end-product, the esti- 
mates of the jet power in the nine elliptical galaxies, is rea- 
sonable. The impact of underestimating the total jet energy 
very nearly cancels out the impact of underestimating the 
timescales involved and the analysis that we have presented 
thus far will only be minimally affected. 

4.3 Reconsidering ADAF models: the effect of 
accretion disk winds 

Finally, we have alluded previously to the possibility that 
the accretion flow may be modified by mass loss in the form 
of winds. Winds from the accretion flo w are likely to b e the 
norm rather than the exception (e.g., lBlandfordll2005h and 
these winds will remove mass, angular momentum and en- 
ergy from the flow. As a consequence, ADAF-like accretion 
flows with winds (advection-dominated inflow-outflow solu- 
tion, ADIOS model, BB99) will have a different dynamical 
structure than ADAFs without winds. For instance given 
the same value of M near the black hole, the angular veloc- 
ity, scale height, total pressure and magnetic field strength 
predicted by the ADAFs with and without winds will be 
different. And since our models for the jet power are linked 
to the accretion flow model through the quantities H, fi' 
and B, it is not inconceivable that jet powers linked to no- 
wind ADAFs could be very different from jet powers linked 
to ADIOS models. 

The ADIOS solution has three parameters in addition 
to those characterizing the no mass loss ADAF model: p w , 
X w and e w . These parameters describe how much mass, en- 
ergy and angular momentum, respectively, the wind removes 
from the accretion flow and specific choice of values for these 
three parameters leads to different types of winds (for more 
information see BB99). To assess the implications for jet 
power if the true underlying model is more correctly de- 
scribed by the ADIOS solution, we begin by fixing the pa- 
rameters common to the two models (the ADIOS and the no- 
wind model) as follows: a — 0.1, 7 = 1.5, and M. = 1O 9 M0 
(we have confirmed that varying these particular variables 
does not affect our conclusions). We then consider various 
ADIOS solutions arising from varying the wind parameters 
across the following range: p w = — 1, X w = 0.1 — 0.75 and 
e w = 0.1 — 0.5. This range of values encompasses all the 
interesting types of winds. 

We found that while some features of the two models 
are very different, the solutions for H , Q' and B at the ra- 
dius of marginally stable orbit in the ADIOS model are, to 
first order, comparable to the corresponding solutions de- 
rived from the no-wind model. In detail, the differences are 
such that the ADIOS jet powers are always smaller than 
the no-wind values by factor of order unity, regardless of 
the combination of values of the wind parameters. Lower 
jet powers in the ADIOS model, in turn, imply both a nar- 
rower range of spins for jet-emitting black holes and that 
these black holes must be spinning faster than those embed- 
ded in the corresponding no mass loss ADAF in order to 



reproduce the correlations derived by A06. In other words, 
if winds are indeed the norm, then the central black holes 
must be spinning even more rapidly than suggested by our 
analysis above. 



4.4 Comparison with previous estimates of jet 
efficiency 

A06 reported an efficiency of conversion of MBondi into jet 
power of Pjet / MsoudiC 2 « 2% for Pj et = 10 43 erg s~ x . In our 
modelling, this corresponds to the case where tBondi = 1 (i.e. 
M ms = MBondi)- Since it is more likely that M ms < M Bo ndi, 
our models suggest that the efficiency of conversion of the 
accreting matter into jet power is considerably higher. For 
instance, rjj ct may reach values as high as « 50% for high 
spin rates. The extraction of spin energy from the holes is 
responsible for this noticeable increase in the jet efficiency. 
We note that the increase in the jet power with black hole 
spin has been ve r ified in numerical simul a tions of jets (e.g., 
De Villiers et al. I 120051 : iMcKinnevI 120051 : lHawlev fc Krolikl 
2006T I. 

The uppper limit we have obtained for r]j ct is more 
than an order of magnitude larger than that reported by 
lArmitage fc Nataraianl (l999). This result is mainly due to 
the fact that we hav e included a Kerr metric shear-driven 
dynamo ijMeierll 19991 ). which enhances the field strength and 
was not included by these a uthors (see q2.1land Appendix ) . 
Based on the cal c ulatio ns of lArmitage fc Nataraianl (|l999h . 
ICao fc Rawlhigsl [|2004r ) have recently asserted that accre- 
tion flow-jet models of the kind we have considered cannot 
account for the jet powers of a sample of low accretion rate 
3CR FR I radio galaxies observed using the Hubble Space 
Telescope (jet powers in the range ~ 10 41 — 10 45 erg s _1 ). As 
a result of much higher efficie ncies, adopting t he critical ac- 
cretion rate M cr it ~ a 2 MEdd IjEsin et al.lll997l ) above which 
the ADAF solution ceases to be valid, and taking j = 0.998 
and a = 0.3, our jet models yield Pjet ~ 10 45-46 erg s _1 , 
which is more than enough t o account for the obser ved jet 
power in the sample used bv lCao fc Rawlingsl (120041 ). 

The adopted model for the jet-formation mechanism 
predicts a linear relation between accretion rate and jet 
power, while the measured best-fit slope of the relation mea- 
sured by A06 is marginably better fit by a non-linear rela- 
tionship. More than likely, this marginal discrepancy is due 
to the small number of systems in the A06 sample. However, 
we note that if tBondi has a slight depend ence on j, as sug- 
gested by recent numerical simulations of lHawlev fc Krolikl 
(2006), then tBondi will also exhibit a weak dependence on 
Pjet and the agreement of the model with the observed cor- 
relation will be improved. We verified that if esondi °c -Pjet 2 
the slope predicted by the model agrees with the observed 
best-fit slope. In the numerical simulations, the dependence 
of Mms on j arises due amplified magnetic fields that act to 
transfer angular momentum from the hole to the accretion 
disk. In addition, it is not inconceivable that the dependence 
eBondi(-Pjct) might also emerge as a consequence of feedback 
effects of the jet on the material fueling the accretion flow. 
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5 CONCLUSIONS 

We have employed two physical models for the black hole- 
accretion flow-jets with features that are broadly consistent 
with the results of numerical simulation to understand the 
empirical correlation between accretion rates and jet pow- 
ers of X-ray luminous elliptical galaxies derived by A06: the 
classical Blandford-Zn ajek model and a hybrid model pro- 
posed by iMeierl (|200ll ). In the BZ model the energy is ex- 
tracted electromagnetically from the spinning black hole; in 
the hybrid model the magnetic field threads both the inner 
and the outer regions of the accretion flow and the energy 
that powers the outflow is extracted from both the rota- 
tion of the disk as well as from the rotation of the black 
hole, via the frame-dragged accretion flow inside the hole's 
ergosphere. We assume that the accretion flow is advection- 
dominated (ADAF) and take into account general relativis- 
tic effects not fully appreciated before in these models. 

In the absense of disk winds and feedback, the model 
suggests that the jet power should be linearly dependent 
on the accretion mass flow rate. The normalisation of this 
relation is dependent on the disk viscosity and the black hole 
spin, with the dispersion around the best-fit correlation of 
A06 being caused by different values of the black hole spin, j, 
and the ratio of the disk mass flow rate to the accretion rate 
onto the black hole, eaondi- We find that the jet efficiency 
(Tjjet = Pjct/MmsC 2 ) can exceed 10% so that the jet may 
carry away an appreciable fraction of the rest mass energy 
of the accreted matterial. 

We compared our jet power models to the jet power esti- 
mates made by A06. Adopting typical values of the viscosity 
parameter a ~ 0.04 — 0.3, the Meondi vs. Pj e t correlation im- 
plies a narrow range of spins j ~ 0.75 — 1 and accretion rates 
M ma m (0.04— l)AfBondi- If we further demand that the mass 
accretion rate be restricted to M ms < Madaf ~ aMBondi, as 
is both likely to be the case and consistent with ADAF model 
expectations, we find that the observed correlations require 
j > 0.9; i.e., the correlations imply rapidly spinning black 
holes in all of the target galaxies. If additionally, the ADAF 
accretion flows also experience mass, energy and angular 
momentum loss via winds as in the ADIOS model proposed 
by BB99, the correlations indicate nearly- maximally spin- 
ning black holes (see i)4.3|l . It is reassuring that semi-analytic 
cosmological simulations of the spin e volution of black holes 
through mergers and gas accretion (|Volonteri et al.l I2005T ) 
and estimates of the radiative efficiencies of global pop- 
ulations of quasars based on Soltan-type arguments (e.g., 
ISoltanl QUI lYu fc Tremamel 12002] ; IWang et al l I2006T ) sug- 
gest that most nearby massive holes are rapidly rotating. 

The relatively small scatter in the correlation between 
^Bondi and Pjct of A06, combined with the strong depen- 
dence of jet power on black hole spin implies large black hole 
spins is probably a general result valid for all elliptical galax- 
ies. Of course the sample considered by A06 is relatively 
small, but if this result holds-up in larger samples, and can 
be shown to also apply to the larger c entral radio galaxies 
of galaxy clusters (IRaffertv et al.l 12006) . our results suggest 
that they have the most powerful jets because they have 
sufficient black hole mass to host relatively high mass flow 
rates while the disk remains in the ADAF state (and hence 
is able to generate strong poloidal magnetic fields near the 
black hole). In this picture, all central ellipticals have the ca- 



pacity to produce powerful jets, and the resulting jet power 
is determined by the structure of the accretion disk. Our 
results reveal a potentially fundamental connection between 
black holes and the formation of the most massive galax- 
ies: the central holes in the cores of galaxy clusters must 
be rapidly rotating, in order to make the radio jets power- 
ful enough to provide an effective feedback mechanism and 
quench the cooling flows, therefore preventing star forma- 
tion and explaining the observed galaxy lum inosity function 
(e.g., ICroton et aLll2006l ; iBower et aLll2006f ) as well as ac- 
counting for the observed X-ray luminosity - temperature 
and X-ray luminosity - Sunyaev-Zeldo vich effect correla- 
tions for the intraclu ster medium (e.g., iBabul et alj 1200*21 ; 
iMcCarthv et alj|2004 ). 
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APPENDIX A: DERIVATION OF THE JET 
POWER 

We list the equations we used to compute the dependence 
of the jet power on a, j and M using the Blandford-Znajek 
model f £|2.2[l and the hybrid model ( i]2.3fl . The jet power 
is given by Equations Q] (BZ model) and [5] (hybrid model) . 
The code that implements the equations described in this 
work and returns the jet power is available at the URL 
http:/ /www. if.ufrgs.br/~rns/jetpower. htm 

The following equations describe the self-similar ADAF 
structure (|Naravan fe Yill 19951 ). where we use the black hole 
mass in solar units (m = M./Mq), accretion rates in Ed- 
dington units (rh = M/M-^dd, M^dd is the Eddington ac- 
cretion rate defined in Q an d radii in Schwarzschild units 
(r = R/{2GM./c 2 )): 

tt' = 7.19 x W 4 c 2 m~ 1 r~ 3/2 s~\ (Al) 

B = 6.55 x 10 8 a" 1/2 (l - P) 1/2 C ; 1/2 cl /4 m' 1/2 rh 1/2 r- 5/4 G, 

(A2) 

H/Rn (2.5c 3 ) 1/2 . (A3) 
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The constants Ci, ci and C3 are denned as 
5 + 2e' 



Cl 



C2 = 



3a 2 
2e'(5 + 2e' 



9'(0L,e') 



nl/2 



9q 2 

2 / / 

. 1 / 5/3-7 
/ V 7-1 



s'(a.e') 



1 + 



18a 2 



1/2 



(5 + 2e') 

The relations among a, 7 and /? are given in t|2.1l The 
angular velocity of the field seen by an outside observer 
at infinity in the Boyer-Lindquist frame is = Q' + w, 
where the angular v elocity of the local metric is given by 
(|Bardeen et al.lll972ft 

_ g$t _ 2aM. 

^ = ~Hm> ~ a?{R + 2M.)+R?' 



(A4) 



using geometrized units [G = c = 1). 

We estimate the field-enhanc ing shear caused by the 
Kerr metric following iMeierl |2001) as g — Q/fi', such that 
the azimuthal component of the field is given by = gB 
(see i|2.ip . The poloidal com ponent is re l ated to the az- 
imuthal component following iLivio et al.l l|l999f) as B p w 
H/R Brf, ~ B,),. 

Lastly, we insert all the quantities defined above into 
Equations[T]and[21 and then evaluate the resulting equations 
at the m arginally stable orb it of the accretion disk 7? ma , 
given by (|Bardeen et allll972l ) 

t = M. {3 + Z 2 - [(3 - Zi)(3 + Zi + 2Z 2 )] 1/2 } , (A5) 



Zi 



1 + (1-J" 

Z 2 



-)Va [(l+i^ + a-j) 1 / 3 ] 



(3/ + Z?) 1 / 2 . 



Taking i? = i? ms in the BZ model corresponds to assume 
that the strength of the magnetic field at the horizon of 
the hole is very similar to the corresponding strength at 
the marginally stable orbit. This is a reasonable assumption 
accor ding to recent numerical simulations of jet for mation 
(e.g., iHirose et alll2004 iMcKinnev fc Gammiell2004h . 



This paper has been typeset from a TjnjX/ I^Tj^X file prepared 
by the author. 



